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Abstract 

A simple trace-coupled Galileon model is shown to admit spherically symmetric static solutions with 
naked spacetime curvature singularities. 



1 



Galileon theories are a class of models for hypothetical scalar fields whose Lagrangians involve 
multilinears of first and second derivatives, but whose nonlinear field equations are still only second 
order. They may be important for the description of large-scale features in astrophysics as well 
as for elementary particle theory 0, 0] ■ Hierarchies of Galileon Lagrangians were first discussed 
mathematically in Qj. The simplest example involves a single scalar field, (f>. This Galileon field 
is usually coupled to all other matter through the trace of the energy- momentum tensor, 

©(matter) 

and is thus gravitation-like by virtue of the similarity between this universal coupling and that of 
the metric g^ u to 0[t™ attcr ^ in general relativity. In fact, some Galileon models have been obtained 
from limits of higher dimensional gravitation theories [6|] . 

In [3] the effects of coupling a Galileon to its own energy-momentum trace were considered, in the 
flat spacetime limit. Here, general relativistic effects are taken into consideration and additional 
features of this same model are explored in curved spacetime [5j. While this investigation was 
in progress, I learned of other work [l| for a related class of models. However, one of the main 
points of that other study is at variance with the results to follow, namely, the model discussed 
here admits solutions with naked singularities when the energy in the scalar field is finite and not 
too large, and for which the effective mass of the system is positive. Thus for the simple model at 
hand there is an open set of physically acceptable scalar field data for which curvature singularities 
are not hidden inside event horizons JljjJnL This would seem to have important implications 
for the cosmic censorship conjecture 0j~12, 15]. It is worthwhile to note that, in gene ral, naked 
singularities have observable consequences that differ from those due to black holes 14]. 

The scalar field part of the action in curved space is 
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This gives a symmetric energy-momentum tensor Q a R for 



— g d 4 x . 
upon variation of the metric. 
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It also gives the field equation for eft upon variation of the scalar field, £ [4>] = 0, where 
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£ [4>] = da 



(4) 
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Since £ [eft] is a total divergence, it easily admits a first integral for static, spherically symmetric 
configurations. Consider only those situations in the following. 

For such configurations the metric in generalized Schwarzschild coordinates is [L 



(ds) 2 = e N ^ (dtf - e L(r) {drf - r 2 {dO) 2 - r 2 sin 2 9 (dip) 2 . 



(6) 



Thus for static, spherically symmetric 0, with covariantly conserved energy- momentum tensor ([5]), 
Einstein's equations reduce to just a pair of coupled lst-order nonlinear equations: 



r 2 9 4 * = e~ L (rV - l) + 1 , 
r 2 9 r r = e~ L (-rN' - l) + 1 . 



(7) 
(8) 
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These are to be combined with the first integral of the <f> field equation in this situation. Defining 

rj (r) = e - L M/ 2 , zu{r)=r, (r) (r) , (9) 

that first integral becomes 

°~ V2 (l + ^) + J(V + -W' (10) 



~^^ = ™\- ■ - , ■ 2 V ■ ,. 

where for asymptotically flat spacetime the constant C is given by lim r _ > . 0O r 2 4>' (r) = C . Then 
upon using 

e t l = @ e e = QV = \zu 2 (1 + \zu 2 ) - r,zu 2 zu' , (11) 

9/ = -\w 2 (1 + \zu 2 ) - i W 3 (N' + I) , (12) 

the remaining steps to follow are clear. 

First, for 0, one can eliminate N' from (jSJ) and (|1U|) to obtain an exact expression for N in 
terms of rj, zu, and C: 

e N/2 = 8C rj-^rzu 3 

rzu (Azu — 2r 2 zu ?J — r 2 zu 5 + 8rr/ + \2wrj 2 + 8rzu 2 r]) 

If the numerator of this last expression vanishes there is an event horizon, otherwise not. When 
n = \rzu^ the denominator of f 1 1 3 [) is positive definite. 

Next, in addition to (|7|) one can now eliminate ./V from either ([8]) or (jlOp to obtain two coupled 
first-order nonlinear equations for r\ and zu. These can be integrated, at least numerically. Or 
they can be used to determine analytically the large and small r behaviors, hence to see if the 
energy and curvature are finite. For example, again for asymptotically flat spacetime, it follows 
that 

e L / 2 ~ i .4- — + - i a\ f2 _ ( \±_ + 1. \ i i r, \ / - or o ) 



r— >oo 



1 + ^ + ^ {6 M 2 -C 2 )^ + 1 -M(5M 2 -2C 2 )^ + 0^ , (14) 

* m ^ l -T-\ M2 V 2 + Y2 M ^ ~ 6M2 ) h + (^) ' (15) 



r— >co 



- ~ ^(l + ^- + lM^)+0(-) . ,10) 



r— >oc t 2 V r 2 r 2 } V 
for constant C and M. 

As of this writing the details of the two remaining first-order ordinary differential equations are 
not pretty, but the equations are numerically tractable. In terms of the variables defined in ([U]), 
Einstein's equation (J7J) becomes 

I (r,zu,rj) r-^-zu + J {r,zu,rj)r^-rj = K (r,w,rj) , (17) 
ar dr 

I (r, zu, rj) = rnzu 2 , J (r, zu, rj) = -2r\ , (18) 
K(r,zu,rj) = \r 2 zu 2 (l + \zu 2 ) + rj 2 - 1 . (19) 

But worse than that, in light of (113p Einstein's equation ([8]) becomes 



F (r,zu,r])r—zu + G (r,zu ,rj)r^-r\ = H (r,zu,n) , (20) 
ar ar 
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F (r, m, rj) = -4rj [2r 3 uj 6 + 3r 3 ro 8 + 16zun + 4rw 4 

+ 16rr] 2 + 48mrj 3 + 48rw 2 n 2 + 12rroV - I2r 2 w 5 rj] , (21) 
G (r, w, rj) = 8nw 2 [2r 2 w 2 + 3r 2 w 4 - 12n 2 + \2rvo 3 n + 4] , (22) 
H (r, w, rj) = w [8r]xu (4rw 3 - 4r] + 2r 2 w 2 r] + 3r 2 w 4 rj + 12rro 3 ?? 2 - I2n 3 ) 

+ (4 + 3r 2 ro 4 + 2r 2 ro 2 + 12?? 2 ) (4-cu - r 2 w 5 - 2r 2 w 3 + 8rvj 2 n + 8rn + 12mi] 2 )] . (23) 

As a representative example with w > 0, ()17p and (|20p were integrated numerically to obtain 
the results shown in Figure 1, for data initialized as vj\ t=1 = 0.5 and r]\ r=1 = 1. Evidently it is 
true that r\ (r) ^ ^ rro3 ( r ) for this case, so e N ^ does not vanish for any r > and there is no 
event horizon. 

However, there is a geometric singularity at r = with divergent scalar curvature: 
lim r ^o r 3 / 2 R = const. Since R = and lim^oro is finite, this divergence in R comes 

from the last term in (|12p . which in turn comes from the second term in A, i.e. the covariant 
d(pd<ftd 2 '4> in (pQ). In fact, it it not difficult to establish analytically for a class of solutions of the 
model, for which the example in Figure 1 is representative, the following limiting behavior holds. 



lim ( e L/2 /^ ) = I , lim ( ^e N/2 ) = n , hmw = p, lim (</>'/ y/f) = pi , (24) 

r — ifl \ / f — \C\ \ J f — ^0 i> — ^0 

the constant C in (Tl0]l : 

2C = 3np 2 /£ . (25) 



r— >0 \ J r— >0 \ J r— >0 r— >0 

where £, n, and p are constants related to the constant C in (I10p : 

,2 



It follows that for solutions in this class, 

lim r 3/2 R = pC/n . (26) 

! — >0 

For the example shown in Figure 1: i ~ 1.5, n ~ 0.086, p ~ 3.3, C ~ 0.94, and pC/n ~ 36. 

For the same ^| r=1 = 1, further numerical results show there are also curvature singularities 
without horizons for smaller ro| r=1 > 0, but event horizons are present for larger scalar fields 
(roughly when ro| r=1 > 2/3). A more precise and complete characterization of the data set 
{ro| r=1 , r/| J1=1 } for which there are naked singularities is in progress, but it is already evident from 
the preceding remarks that the set has nonzero measure. 

The energy contained in only the scalar field in the curved spacetime is given by 



£ G aiiieon= / H (r) dr = e s %{e s )ds, (27) 

j0 J -oc 

U (r) = 4i:r 2 e L l 2 e N l 2 Q t t = 2^e 2s e L l 2 e N ' 2 w 2 (s) (l + \w 2 (s)) - 4^e s e N ' 2 w 2 (s) £w (s) . (28) 



For the above numerical example, the integrand e s % (e s ) is shown in Figure 2. Evidently, -Ecaiiieon 
is finite in this case. It is also clear from the Figures that the Galileon field has significant effects on 
the geometry in the vicinity of the peak of its radial energy density. There the metric coefficients 
are greatly distorted from the familiar Schwarzschild values, and as a consequence, the horizon is 
eliminated. 

It remains to investigate the stability of the static solutions described above, and to consider 
the dynamical evolution of generic Galileon and other matter field initial data, along the lines of 
0] , to determine under what physical conditions the naked singularities discussed here are actually 
formed. 

/ thank D Fairlie for many discussions about Galileons in general, and this work in particular. I 
also thank S K Rama, N Rinaldi, R Wald, and K S Virbhadra for discussions of naked singularities. 
This research was supported by a University of Miami Cooper Fellowship, and by NSF Awards 
PHY-0855386 and PHY-12U521. 
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Fig. 1: For initial values w(s)\ s=0 = 0.5 and r](s)\ s=0 = 1.0, d<p/dr = w/r] is shown in red, 
e = 1/rj 2 in green, and e N in blue, where r = e s . For comparison, Schwarzschild e L and e N are 
also shown as resp. green and blue dashed curves for the same M ~ 0.21 
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